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ABSTRACT 


The  thrust  of  this  paper  is  to  address  the  question  of  what 
dynamical  properties  of  a viscoelastic  liquid  lubricant  are  instru- 
mental in  determining  the  traction  in  a heavily  loaded  rolling  contact. 

In  examining  this  question  we  first  explore  the  dynamical  changes  that 
the  large  normal  stresses  in  the  contact  zone  produce  in  the  lubricant 
film.  The  effect  of  these  structural  changes  on  the  response  of  the 
liquid  to  the  imposed  shear  is  then  evaluated;  from  this,  the  shear 
stress  profile  in  the  contact — which  determines  the  traction — can  be 
computed. 

The  principal  features  of  the  behavior  can  be  summarized  as 
follows.  As  a fluid  element  moves  through  the  contact  region,  the  large 
and  rapidly  changing  pressure  to  which  it  is  subjected  produces  alter- 
ations of  the  liquid  structure  that  reduce  the  molecular  mobility. 
Depending  upon  the  detailed  nature  of  the  structural  relaxation  dynamics, 
this  leads  to  a glassy  or  solid-like  behavior  over  at  least  part  of  the 
zone.  As  a result  the  shear  response  characteristics  (which  depend  on 
the  instantaneous  structure)  vary  as  the  element  moves  through  the  con- 
tact. The  effect  is  that  the  stress  profile  is  intermediate  between 
those  typical  of  low  viscosity  liquids  and  amorphous  solids.  For  a given 
set  of  external  variables  (load,  rolling  speed)  the  initial  slope  of  the 
traction  curves  is  found  to  be  governed  mainly  by  the  forms  of  the  shear 
and  structural  relaxation  functions  and  the  shear  rigidity  modulus. 

Other  properties  such  as  the  inlet  viscosity  and  the  pressure-viscosity 
coefficient  are  found  to  play  somewhat  less  important  roles. 


I . INTRODUCTION 


Traction  is  the  transfer  of  power  from  one  rolling  element 
(the  driver)  to  a second  such  element  (the  slave)  through  the 
shearing  of  the  lubricant  film  that  separates  the  two  pieces.  In 
elastohydrodynamic  lubrication  (EHDL) the  lubricant  contact  is 
heavily  loaded  so  that  the  elastic  deformation  of  the  rolling  mem- 
bers is  not  insignificant  compared  to  the  thickness  of  the  lubri- 
cant film.  In  recent  years,  EHDL  theory  has  progressed  to  a ra- 
ther sophisticated  stage  but  as  yet  is  incapable  of  yielding  quan- 
titatively meaningful  results  for  the  traction  between  driver  and 
slave.  The  principal  difficulty  arises  when  one  attempts  to  de- 
scribe the  lubricant  in  the  contact  zone. 

In  this  paper  we  address  directly  the  two-fold  question: 

How  does  a liquid  lubricant  behave  in  a heavily  loaded  rolling  con- 
tact and  what  are  the  fundamental  properties  of  a liquid  that  are 
instrumental  in  determining  the  traction  in  such  circumstances. 

For  definiteness  we  examine  the  situation  exemplified  by  the  "twin- 
disk"  apparatus  commonly  used  in  measuring  the  traction  in  a line 
contact.  Fig.  1 gives  a schematic  representation  of  the  geometry. 

The  quantity  of  interest  is  the  traction  coefficient 

CT  = Ft/Fn 

where  Ft  is  the  tangential  force  in  the  contact  (the  tractive  force) 
and  Fn  is  the  normal  force  or  load  exerted  perpendicular  to  the  con- 
tact. 

In  the  next  section  of  this  paper  we  analyze  the  behavior  of 
the  liquid  in  the  contact  zone  and  establish  a theoretical  frame- 
work for  calculating  the  traction  in  terms  of  basic  liquid  properties. 


lubricant  parameters  and  give  a brief  discussion  of  the  signifi- 


cance of  these  results. 


3. 


II  THEORETICAL  FORMULATION 

Qualitatively,  the  picture  of  lubricant  behavior  in  an  EHDL 
contact  is  rather  straightforward.  Under  the  high  pressures  ex- 
perienced in  the  contact  zone,  the  mobility  of  the  molecules  in 
the  liquid  is  severely  restricted.  Consequently  for  much  of  the 
time  required  for  a fluid  element  to  traverse  this  zone,  the  li- 
quid structure  is  dynamically  arrested  or  "frozen  in;  as  a result 
the  lubricant  exhibits  glassy  or  amorphous  solid-like  character- 
istics. This  solid  pad  separating  the  rolling  elements  transmits 
the  tangential  traction  force  through  the  essentially  elastic 
shearing  stresses  developed  in  it.  The  theoretical  formulation 
developed  in  this  section  translates  these  ideas  into  a mathemati- 
cal framework  by  which  the  stress  profile  in  the  lubricant  film 
and  thus  the  tractive  force  can  be  calculated. 

1.  Preliminaries 

The  contact  zone  geometry  and  the  coordinate  axes  are  shown 
in  Fig.  2.  The  disks  rolling  with  velocities  U^  and  U2  (Ui  > U2) 
have  deformed  elastically  to  form  an  EHDL  contact  zone  of  (approxi- 
mately) constant  thickness  2h  and  width  2b.  The  pressure  distri- 
bution in  the  zone  is  assumed  to  be  uniform  in  the  y-direction 
and  to  vary  with  x as  an  elliptical  Hertzian  profile 

p(x)  * PHzri-(x/b)2]l/2  (|x|  4 b),  (2) 

where  PHZ  is  the  maximum  Hertzian  pressure.  The  half-width  of  the 

(2) 


zone  is  given  by 
b = 4RPHZ/E, 


(3) 


where  R = RiR2/CRi  + R2)is  an  effective  radius  for  the  disk  pair 
of  radii  R^  and  R2.  Here  E is  an  effective  elastic  modulus  for  the 
disks  related  to  the  Young's  modulus  Eq  by  E = Eg/(l-v2)  where  v 
is  Poisson's  ratio.  The  normal  force  per  unit  (z  - direction) 
length  is  thus  easily  calculated: 

fn  * dx  P(x*  = irfa  pHz/2-  (4) 

Understanding  and  calculating  traction  basically  means  ex- 
plaining quantitatively  the  so-called  "traction  curves" — plots  of 
traction  coefficient  versus  slide-roll  ratio  AU/U — such  as  those 
shown  in  Fig.  3.  The  sliding  (or  slip)  speed  AU  = Uj,  - U2  and  the 
mean  rolling  speed  U = (U]_  + U2)/2.  In  this  paper  we  shall  restrict 
our  consideration  to  the  linear  portion  (the  low-slip  region)  of 
these  curves . 

The  analysis  is  simplified  further  by  the  following  approxi- 
mations: (1)  All  the  quantities  are  independent  of  z,  the  co-ordi- 
nate parallel  to  the  axes  of  rotation  of  the  disks;  (2)  there  are 
no  temperature  gradients  within  the  lubricant  in  the  contact  zone; 

(3)  the  only  velocity  gradients  in  the  liquid  that  must  be  con- 
sidered are  those  in  the  y-direction;  and  (4) all  effects  arising 
from  the  curvature  of  the  disks  or  irregularities  in  the  film 
thickness  can  be  ignored.  The  problem  is  thus  one  of  computing 
the  shear  stress  distribution  within  a liquid  of  rectangular  cross- 
section  (2b  x 2h)  subjected  to  a Hertzian  pressure  distribution  and 
a constant  shearing  rate 

e = 3u/3y  = U/2h,  (5) 

where  u is  the  x-component  of  the  particle  velocity. 


5. 


2.  Tractive  Force 

The  tractive  force  per  unit  (z-direction)  length  ft  trans- 
ferred across  the  liquid  film  is  given  by  the  integral  of  the 


shear  stress  profile  CyXCx)  over  the  EHDL  contact  zone: 


,+b 


f +-  — / 1-  dX  tjyy  (.X)  • 
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Rather  than  examining  the  entire  lubricant  film  at  a fixed  instant, 
as  implied  by  Eq.  (6) , it  is  somewhat  more  convenient  to  consider 
a fluid  element  of  dimensions  dx  *2h  travelling  through  the  con- 
tact (from  x « -b  to  x = +b)  at  a speed  U.  Formally  this  amounts 
to  using  the  transformation  x -*■  -b  + Ut  (where  x locates  the  posi- 


tion of  the  element  at  time  t)  in  Eq.  (6) : 
t0 

f+.  = U /n  dt  a (t ) , 


(7) 


where  a (t)  = CyX  (-b  +Ut)  is  the  shear  stress  on  the  element  at 
time  t and 

t0  = 2b/U  (8) 

is  the  transit  time  of  the  fluid  element  through  the  contact  re- 
gion. The  problem  is  thus  to  construct  the  form  of  the  time  de- 
pendent stress  a(t). 

3.  Linear  Viscoelastic  Theory 

Because  we  are  dealing  only  with  small  shearing  rates,  we 
shall  employ  linear  viscoelasticity  to  describe  the  shear  stress- 
strain  relationship.  The  most  general  form  allowable  is ^ 


cx(t)  = / 0 dt7  $s  (t-t' ) e (t') . 


(9) 


Equation  (9)  gives  the  shear  stress  at  time  t as  a summation  (or 
integral)  of  contributions  arising  from  shear  strains  e(t')dt' 
imposed  on  the  fluid  at  times  0 ^ t' 4 t (it  is  assumed  that  no 


' ^ * » V' ak  . . 
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strains  are  present  prior  to  time  t7  = 0 and,  because  of  causality 
no  strains  imposed  after  ty  = t can  contribute  to  the  stress  at 
time  t)  . The  shear  stress  relaxation  function  SgCt-t7)  gives  the 
relative  weightings  at  t of  the  various  contributions;  that  is, 

4>s(t  -t-'JeCt')  dt7  is  the  stress  at  time  t resulting  from  a strain 
e (t7 ) dt7  imposed  at  time  tt 

This  can  be  clarified  a bit  more  by  examining  several  familiar 
cases.  A simple  Newtonian  viscous  fluid  is  characterized  by  the 
fact  that  there  is  an  instantaneous  proportionality  between  the 
stress  and  the  shear  rate.  This  is  expressed  by  writing  $s  as 
a constant  n multiplied  by  6 (t  - t'),  the  Dirac  delta-function; 
that  is 

*s(t  - t')  = n <$  (t-t*)  . (10) 

If  we  use  this  form  in  Eq.  (9),  the  result  is 

a (t)  = ne (t) , (11) 

the  familiar  Newtonian  form  with  n being  the  shear  viscosity. 

The  case  of  a linear  elastic  solid  in  which  there  is  a con- 
stant proportionality  between  the  stress  and  the  strain  is  de- 
scribed by 

$s(t-t/)  = Gm  (t  » t') . (12) 

In  this  case  Eq.  (9)  leads  to 

t , . / 

a(t)  = G,,,  fQ  dt  e(t)  = Gwe(t).  (13) 

Here  £(t)  is  the  shear  strain  and  thus  g»  is  shear  rigidity 
modulus . 


For  a viscoelastic  material  the  form  of  $s(t-t' ) is  inter- 
mediate between  the  instantaneous  decay  ($s  = 0 for  t>t') 


7. 


appropriate  for  the  Newtonian  fluid  and  the  constancy  exhibited 
for  an  elastic  solid.  The  form  must  be  such  that  the  response 
of  the  system,  i.e.,  the  stress,  to  an  applied  strain  is  (a)  purely 
elastic  immediately  following  the  strain  application  (t-ty -*•()),  and 
(b)  purely  viscous  after  a long  time  has  passed  (t-t+  ») . To 
accomodate  the  first  of  these  requirements  is  is  sufficient  to  de- 
mand that 

♦a  CO)  = 0.,  U4> 

0 

while  to  take  care  of  the  second,  it  is  necessary  that 

$s  (t-t-*- °°) -*■  0.  Cl 5 ) 

The  vanishing  of  <J>S  at  large  values  of  its  argument  must  be  suf- 

t f t 

ficiently  rapid  that  / g dt*  $s  (t)  is  finite.  That  this  must  be  so 

can  be  seen  from  the  following  considerations.  Suppose  that  a 

shear  rate  e is  imposed  on  the  system  at  time  zero  and  maintained 

constant  thereafter.  From  Eq.  (9) 

. t . t 

a(t)  = e /0  dt'  4>s  (t-t' ) = e /q  dt"  $s(t")f  (16) 

where  the  last  equality  follows  from  the  substitution  t = t-t  . 
Taking  the  t-*»  limit  and  comparing  the  result  with  Eq.  (11)  yields 
/Q  dt  *s(t)  -n?  U7) 

the  integral  is  indeed  finite  and  is  in  fact  just  the  shear  vis- 
cosity of  the  liquid. 

It  is  convenient  to  construct  a normalized  stress  relaxation 
function  <j>g  (t)  by  the  definition 

*SCt)  = G^s  (t)  . (18) 


Clearly  this  implies 
4>s  to)  = 1 
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(19) 

and 

/0  dt  <t>s  (t)  = n/Goo  H Tg  (20) 

Equation  (.20)  serves  as  the  definition  of  the  shear  relaxation 
time  Tg.  Recalling  the  normalization  condition  (19)  , it  is 
apparent  that  xs  provides  a rough  characterization  of  the  time 
scale  for  $s  to  decay  by  some  significant  fraction  of  its  original 
value. 

Observe  also  that  the  definition  of  ts-'-in  Eq.  (20)  is  the 
equivalent  of  the  familiar  Maxwell  relaxation  time  even  though 
the  definition  is  independent  of  any  specific  viscoelastic  model. 
The  usual  Maxwell  single  relaxation  equations  are  obtained  for 
the  specific  exponential  decay  form  of  s : 

<PS  (t)  = exp  (-tG»/n)  • (21) 

In  passing  we  also  observe  that  the  commonly  employed  fre- 
quency dependent  shear  viscosity  n(w)  and  modulus  G(oj)  are  ob- 
tained from  the  real  and  imaginary  parts  of  the  one-sided  Fourier 
transform  of  $s(t).  That  is, 

n(cj)  - iG (oj) / gj  = fZ  dt  e $s  (t)  . (22) 

With  this  as  background  we  can  now  examine  Eq.  (9)  under  the 
conditions  of  interest  in  this  problem,  namely  when  e (t)  is  a con- 
stant as  indicated  in  Eq.  (5) . In  this  case 

a(t)  = eG«o/|dtn4.s(t")  (23) 

where  we  have  made  use  of  Eq.  (18)  and  have  let  t"  = t-t'.  The 
quantity  GaB/£dt"<Js  (t")  is  clearly  seen  to  be  the  ratio  of  the 


-r,-TOv-r 


instantaneous  stress  to  the  shear  rate;  as  such,  we  can  think  of 

* 

it  as  a time  dependent  effective  viscosity.  It  is  to  be  empha- 
sized that  this  effective  viscosity, 

y(t)  = dt"  <t>s  (t")  , C24) 

is  meaningful  only  under  the  constant  shear  rate  conditions  for 
which  it  is  defined,  and,  moreover,  that  it  approaches  the  actual 
,,true,,viscosity  of  the  liquid  only  in  the  limit  t -*•  <»  as  is  evident 
from  a comparison  with  Eq.  (17) . The  formula  for  ft  can  be  written 
as 

t 

ft  = (Ut0)  [^-  f ° dt  y (t ) ] c|£)  . (25) 

The  factor  UtQ  = 2b  is  just  the  width  of  the  contact  zone;  the 
quantity  in  brackets  can  be  defined  as  an  average  effective  vis- 
cosity 

«:-/  O dt  y(t)  . (26)' 

fco  ° 

In  terms  of  this  and  with  the  help  of  Eq. (4) , the  traction  coef- 
ficient can  be  written  as 

CT  = (2b/irhPHZ)  [Gm(— JL— : ) ] (AU/U).  (27) 

^^00 ’"O 

Apart  from  the  external  parameters  b,  h and  PHZ,  the  slope  of  the 
traction  curves  (Fig. 3) is  proportional  to  (but  not  equal  to) the 
shear  modulus  G^.  The  coefficient  of  proportionality  is  the  dynamic 
factor  iT/JjG^tQ;  as  we  shall  see,  this  factor  never  exceeds  unity. 
Notice  that  since  b * PHZ  the  dependence  of  CT  on  Hertzian  pres- 
sure arises  solely  from  the  dependence  of  Ga  and  y on  the  pressure. 

At  this  point  the  calculation  of  the  tractive  force  would  be 
essentially  completed  if  it  were  possible  to  ignore  the  effects  of 


t 


10. 


the  pressure  pulse  to  which  the  liquid  is  subjected  as  it  passes 
through  the  contact.  One  would  use  the  {presumably  known)  value  of 
Go,  an<^  the  functional  form  of  <{>s(t),  determined  either  experimen- 
tally or  from  some  molecular  model,  and  compute  the  stress  via 
Eq. (23) . This  being  done,  the  tractive  force  and  the  traction 
coefficient  follow  immediately. 

It  is  instructive  to  carry  out  this  program  for  a simple  situa- 
tion ignoring  the  pressure  induced  effects  as  some  useful  limits 
can  be  obtained.  Let  us  assume  that  the  Maxwell  form  of  <t>s(t) 
given  in  Eq.  (21)  is  appropriate.  The  effective  viscosity  is  easily 
calculated  as 

UM(t)  = Tg  [1-exp  (-t/Ts)  ] (28) 

where  we  have  appended  the  subscript  "M"  to  denote  the  use  of  the 
Maxwell  stress  relaxation  function.  Multiplying  this  by  £ =U/2h 
gives  the  tractive  force: 

ft  = (UtoHG^  T s [ 1“  ^ (l-d"to/Ts)]}  (|g).  (29) 

Notice  that  UtQ  = 2b  so  that  the  only  variation  with  rolling  speed 

enters  through  the  dependence  of  p-,  the  quantity  in  curly  brackets, 

on  tQ.  At  low  viscosities  or  low  rolling  speeds  such  that  xs<<t0, 

the  exponential  can  be  ignored  and 

ft  = (2b)  [n(l-Ts/t0)J  (§g).  (x2  » 1)  (30) 

* s 

t 

In  this  limit  the  lubricant  behaves  as  a simple  viscous  liquid.  The 


4.  Structural  Relaxation  in  EHDL 

It  is  clear  from  Eq.  (27)  that  the  crucial  quantity  in  deter- 
mining the  traction  curve  slopes  is  the  factor  v/hG^tQ.  This 
factor  varies  from  2xs/t0,  which  is  essentially  zero,  in  the  low 
viscosity  (ts  <<  tQ)  limit  to  unity  in  the  high  viscosity  C*s  >>  tQ) 
limit.  Because  of  the  large  pressure  pulse  to  which  the  liquid  is 
subjected  in  its  passage  through  the  EHDL  contact,  the  relaxation 
time  varies  over  many  orders  of  magnitude  ranging  from  xs  <<  tQ 
at  the  inlet  to  (and  exit  from)  the  contact  to  ts  >>  tQ  near  the 
center.  It  is  necessary  to  account  for  this  variation  when  calcu- 
lating the  stress  and/or  effective  viscosity  from  Eq.  (23)  and  (24) . 
To  do  this  we  recognize  that  the  relaxation  time  xs  is  determined 
not  only  by  the  instantaneous  pressure  and  temperature  but  also  by 
the  local  structure  of  the  fluid.  It  is  therefore  necessary  to  de- 
vise a theoretical  fomulation  in  terms  of  which  the  structural  dy- 
namics can  be  computed. 

A.  Description  of  the  Dynamical  Structure.  The  thermodynamic 
state  of  a system  in  equilibrium  is  determined  when  the  pressure  P 
and  temperature  T are  specified  (in  a multi-component  system  one 
must  also  give  the  concentrations  of  the  various  species) . We  sup- 
pose that  it  is  possible  to  define  a parameter — for  which  we  use 
the  symbol  C — that  describes  the  local  structure  of  the  material. 

For  a system  in  equilibrium  5=5=  S(P/T).  In  what  follows  we 
shall  be  concerned  only  with  isothermal  processes  and  consequently 
we  shall  suppress  the  explicit  dependence  of  £ on  temperature; 
e.g.,  we  write  £ = E (P)  . 


13. 

If  a change  in  the  pressure  AP  is  imposed  on  the  system  at 
some  time,  say  at  t=t]_,  the  structure  does  not  adjust  instantaneously 
to  the  new  conditions  and  thus  immediately  following  the  impressed 
change  £(t)  E(P+AP);  only  as  t-* 00  can  we  expect  £ (t)  to  reach  its 
new  equilibrium  value.  Practically  speaking,  of  course,  an  infinite 
time  interval  is  not  required  for  structural  equilibration.  The  time 
scale  over  which  the  structure  relaxes  to  its  new  equilibrium  value 
is  described  by  the  structural  relaxation  time . t,  and  the  dy- 
namical  evolution  of  £ is  characterized  by  the  structural  relaxation  — i 
function  ( t/x ) . 

Specifically  consider  the  following  situation.  A system  in  equi- 
librium for  t<t]_  at  pressure  PQ,  so  that 

£(t<tx)  = S(P0)  s £o,  (38) 

is  subjected  to  a pressure  increase  AP  at  t=t]_.  When  t-*- 00  the  system 
will  equilibrate  at  P0+AP,  the  new  equilibrium  value  being 
£(t--  ) = 5 (PQ  + AP)  = £„  . (39) 

The  instantaneous  value  of  £ Is  given  by  an  equation  of  the  form 
£ (t)  =£0+y (£»“£0)  + (1-Y)  (£«-CoHl-<H(t-t;L)/Tl},  (40) 

with  0$Y$1.  Here  the  term  Y(£a,-£0)  gives  that  part  of  the  change 
in  £ that  occurs  promptly  with  the  pressure  step.  The  quantity 
(1-Y)  (£oo~£o)  is  therefore  the  magnitude  of  the  time-dependent  or 
relaxational  change  in  £.  The  factor  l-<>»  [ (t-t^)/T)  gives  the  spe- 
cific time  dependent  approach  to  equilibrium.  As  with  <j>s(cf.  Eqs. 

(15)  and  (19)],  $ is  normalized  so  that  <(>(0)  = 1 and  <J>  (°°)  = 0.  The 


► 


behavior  described  by  Eq.  (40)  is  illustrated  in  Fig.  4.  In  this 
example  we  have  presumed  that  the  pressure  increase  is  sufficiently 
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small  that  the  response  can  be  treated  linearly;  specifically  this 
means  that  x can  be  taken  as  constant.  Regarding  x as  a function 
of  pressure,  we  are  assuming  that  is  sufficiently  small  that 

x (P0+AP)  = x (PQ)  + C3t/3P)p  AP  +...  = x (PQ)  , (41) 

o 

and  thus  one  need  not  worry  about  computing  the  change  in  x oc- 
casioned by  the  pressure  change. 

If  the  imposed  pressure  changes  are  large,  such  as  is  the  case 
in  EHDL  contacts,  then  this  linearization  is  not  possible  and  a more 
detailed  analysis,  which  includes  the  changing  of  the  structural 
relaxation  dynamics "themselves  in  response  to  the  pressure  changes, 
is  required.  The  result  of  this  more  complete  treatment  is  the 
following  generalization  of  Eq.  (40) : 

t 

C (t)  = CQ  + .Y/t  dt'  (df/dt') 

t t 

+ (l-Y)/t.,  dt'  (dC/df)  (1-<J>[A  dt " / x (t")  ] } • (42) 

A heuristic  derivation  of  this  result  is  presented  below. 

Suppose  that  the  system  (in  equilibrium  at  PQ  for  t<t]_)  is 
subjected  to  an  externally  impressed  pressure  variation 

P(t)  = PQ  + P(t-tx)  (43) 

for  t£t]_.  For  any  time  t>t^  we  may  analyze  the  behavior  in  an 
approximate  fashion  by  dividing  the  time  tj_  to  t into  N intervals 
by  a set  of  "cuts"  at  t2 / t3. . .t^. . . tjj.  The  duration  of  any  of 
these  intervals,  e.g.  the  k— , is  denoted  by 

Atfc  = tfc+i  - t^. 


(44) 
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We  also  replace  the  continuous  pressure  variation  pCt-t^  by  a 
series  of  pressure  steps  APj.  imposed  at  t^,  AP2  imposed  at  t2*.. 

The  amplitude  of  the  pressure  step  imposed  at  the  start  of  a given 
interval  is  such  as  to  make  the  pressure  during  that  interval  have 
the  value  that  would  be  given  by  Eq. C43)  for  the  midpoint  of  the 
interval.  Thus  at  time  tk,  the  imposed  pressure  step  is 

APk  = Pl*Ctk+1  + tk)  -t-Ll-ptUtjc  + t^-tj  C45) 

so  that  the  pressure  from  tk  until  tk+]_  -is*'  Po+p  (tk+tk+1)  -tl^  * 

The  approximation  is  illustrated  in  Fig.  5.  Note  that  in  the  limit 
N+®  (Atk->-0)  the  approximation  becomes  exact.  Each  interval  is 
chosen  such  that  the  response  of  the  system  to  any  of  the  pressure 
steps  can  be  treated  linearly,  i.e.  using  Eq.  (40) . Thus  at  time  t ^ 
(the  beginning  of  the  first  interval)  a pressure  P1  = p[4(t2~ti)] 
is  applied;  the  structure  changes  in  response  to  this  so  that  at  t2 
(the  end  of  the  first  interval) 

S(t2)  = C0+YAC1+  (1-Y)  A5i{l-«[(t2-ti)/f1]},  (46) 

where 

ACi  = (3?/3t)t  Ati  = OC/3P)p0  (dp/dt)  tl  Atx,  (47) 

and  t ^ is  the  structural  relaxation  time  for  the  thermodynamic 
and  structural  conditions  appropriate  to  the  tj_  to  t2  interval  (the 
nature  of  the  dependence  of  t on  the  "thermodynamic  and  structural 
conditions"  is  elaborated  below) . At  the  beginning  of  the  second 
interval  (time  t2)  a second  pressure  step 
p2  = Pt*j(t3+t2)-ti]  - p[h(t2~t1)] 

is  applied  and  we  must  now  account  for  (1)  the  continuing  evolution 
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of  the  structure  to  the  first  pressure  step  APj_  except  that  during 
the  time  t2  to  t3,  the  relaxation  time  assumes  a new  value  T2  appro- 
priate for  the  conditions  during  this  second  interval;  and  (2) the 
response  of  the  system  to  the  pressure  change  AP2.  Thus  at  time  t3, 
the  time  dependent  part  of  the  response  to  AP]_  will  be  given  in 
terms  of 

♦ [ (t2~ti)/xi  + Ct3-t2)/T2] 
and  the  response  to  AP2  in  terms  of 

♦ Ct3-t2)/  T 2 ] . 

The  value  of  5 at  t^  is  therefore 

C(t3)  = Co  +YA?!  + (l-Y)AC1{1_4>[Ct2-ti)/T1  + C t3— 1 2/ x 2 ] } 

-hYAC2  + (1-  )AC2{l-<M(t3-t2)/T2]}  (48) 

where  A?2  is  defined  analogously  to  A?i  in  Eq. (47) . 

The  form  of  the  general  result  now  becomes  clear.  At  time  t 
(the  end  of  the  N—  interval)  we  must  account  for  the  response  to  N 
pressure  steps  APi>  AP2...APjj  imposed  at  times  t^/t2...tN  respectively 
The  relaxation  function  describing  the  dynamical  part  of  the  response 
to  the  k—  of  these  pressure  changes  is 

^ f (tk+l“tk^ ^Tk  + (tk+2_tk+l^ /Tk+1  + ,**  + (t“tjj^/xN^ 
where  xk  is  the  structural  relaxation  time  for  the  thermodynamic 

and  structural  conditions  of  the  k—  interval.  Therefore  we  get  the 
result 

N N N 

C(t)=C0+Yk£1ACk+Cl-Y)k|1ACk{l-(J.[j|Jc(tj+1-tj.)/Tj]  . (49) 

N 

Taking  the  limit  at ^0  (Zj=1At  j=t-t]_)  and  letting  Atj-^dt"  , 

t j-*T  (t" ) and  A?^-*- (3C/3t ' ) dt ' allows  us  to  make  the  replacements 


N t 
Z f 
k=l  ti 
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N t 

Z f 
j=k  f 

which  immediately  gives  Eq.  (42)  . An  equivalent  and  often  more 
useful  form  of  this  equation  is  obtained  if  the  second  equality 
in  Eq.  (47)  is  used  to  rewrite  the  A^.  Here 
5 (t)  = q + Y dt' or/3P)  Cdp/dt') 

+ (l-Y)/t1  dt’ OT/3P)  Cdp/dt')  {1-<H/^  dt"/T(t")  ] }.  (50) 

In  order  to  apply  Eq.  (42)  or  (50)  , it  is  necessary  to 
specify  the  nature  of  the  structural  ordering  parameter  the  de- 
pendence of  t on  P and  5 (and  thus  on  time)  and  the  form  of  the 

. 

structural  relaxation  function  ^(t/r).  In  addition  the  interplay 
between  the  structural  and  shear  relaxation  processes  must  be  es- 
tablished. We  shall  now  take  up  these  questions. 

B.  The  Fictive  Pressure.  Under  the  isothermal  circumstances 
under  consideration  here  a convenient  choice  for  the  ordering  para- 
meter describing  the  structure  is  the  fictive  pressure.  A material 
is  said  to  have  a fictive  pressure  Pf  if  the  structure  of  the  ma- 
terial is  the  structure  that  would  obtain  in  equilibrium  at  a 

I 

pressure  P = Pf.  This  definition  is  analogous  to  the  "fictive 

4 c 

temperature"  concept  introduced  and  explicated  by  Tool  and  Ritland 
to  describe  the  structural  changes  in  glasses  as  they  are  cooled 
from  the  melt  through  their  transformation  regions. 


In  the  present  context  the  dynamical  structure  is  described 
at  any  instant  by  a time  varying  fictive  pressure  Pf  (t)  . The  spe- 
cification of  a non-equilibrium  state  is  assumed  complete  when 
Pf(t)  and  the  instantaneous  pressure  P (t)  are  given.  The  dynamical 
behavior  of  Pf(t}  as  it  evolves  toward  its  equilibrium  value  P(t) 
(obviously,  from  the  way  it  is  defined,  the  equilibrium  fictive 
pressure  is  just  Pf  = P) . If  a small  pressure  increase  AP  is  imposed 
at  time  t^  on  a system  previously  in  equilibrium  at  pressure  PQ, 
the  time  variation  of  Pf  (t)  is  given  by  'an  equation  of  the  form  of 
Eq.  (40) . Specifically,  we  take: 

5 (t)  = Pf  (t) . The  structural  ordering  parameter 
is  the  fictive  pressure. 

£0  = P0  and  5^=  PQ  + AP.  The  equilibrium  value 
of  the  fictive  pressure  is  the  actual  pressure. 

y = 0.  Only  relaxational  changes  in  structure 
are  characterized  in  terms  of  Pf;  any  instanta- 
neous changes  are  presumed  to  be  given  explicitly 
in  terms  of  the  change  in  P. 

Thus,  Eq.  (4)  becomes 

Pf(t)  = P0  + AP{1-<H  (t-t^/Tj  }.  (51) 

The  response  to  an  arbitrary  pressure  variation  is  given  by  the 

generalization  of  the  form  of  Eq.  (50) : 

t t 

Pf(t)  = Pg  + /tl  dt ' (dp/dt ' ) {l-<ji  [/fc  idt"/T  (t" ) J } (52) 

where  we  have  used  the  fact  that  3Pf/3p  = 1 to  simplify  the 
integrand. 
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In  order  to  evaluate  this,  the  behavior  of  r with  time  must 

6 7 

be  known.  In  the  spirit  of  Narayanaswamy  and  Gardon  ' twho  deal 
with  isobaric  temperature  changes)  we  can  deal  with  this  phenome- 
nologically by  recognizing  that  t depends  on  both  P and  Pf.  If  the 
structural  relaxation  time  measured  as  an  equilibrium  material  pro- 
perty, varies  with  pressure  as 

Teq  = f(P),  (53) 

then  we  assume  that,  for  a given  non-equilibrium  state  (characterized 
by  P and  Pf) , t is  given  by 

t = f[xP  + (l-x)Pf]  (54) 

where  x (0<x<l)  is  a parameter  characterizing  the  relative  depen- 
dence of  the  relaxation  time  on  pressure  and  fictive  pressure. 

We  designate  it  as  the  "structure  coefficient."  Now  since  both  P 
and  Pf  are  time  varying  we  use  this  in  Eq.  (52)  as 

x(t")  = f[xP(t")  + (l-x)Pf  (t")]  . (55) 

Observe  that  in  order  to  compute  Pf(t)  from  Eq.  (52) , it  is  necessary 
that  T(t") — and  thus  Pf(t") — be  known  only  for  t^<t"5t.  One  can 
thus  proceed  in  a stepwise  fashion  to  calculate  P^(t)  for  all  t>t^ 
since,  presumably,  x(t^)  is  known  as  an  equilibrium  property. 

This  provides  us  with  a clue  for  computing  the  shear  stress 
variation  a(t)  or  equivalently  the  effective  viscosity  u (t)  in  the 
presence  of  the  Hertzian  pressure  pulse.  We  examine  Eq.  (23)  or 
Eq. (24)  and  regard  the  shear  relaxation  function  $s(t")  as  a func- 
tion of  the  reduced  time  variable  t"/xg.  We  then  assume  that  xs 
varies  in  time  through  the  contact  zone  as  a result  of  the  changing 
pressure  and  structure  in  a manner  analogous  to  Eq.  (55)  for  the 
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strucutral  relaxation  time;  viz., 

TsCt")  = g[xP(t")  + U>xl  Pf(t")J.  (56) 

where  the  function  g(P)  gives  the  pressure  variation  of  the  shear 

relaxation  time  measured  as  an  equilibrium  property.  Using  this 

the  effective  viscosity,  for  example,  would  be  given  by 

t t 

U (t)  = Ga,  ftl  df  *aUt'  dt"/Ta  (t”)  ] (57) 

where  the  generalization  of  the  argument  of  <j>s  comes  about  in  es- 
sentially the  same  way  that  the  similar  generalization  of  the  ar- 
gument of  <t>  does  in  Eq.  (50)  . Note  that  in  writing  Eg.  (57)  we 
have  assumed  that  the  variation  of  G„  in  the  contact  zone  can  be 
ingored.  This  is  not  a necessary  restriction;  indeed  by  characteri- 
zing the  dependence  of  Goo  on  pressure  and  structure,  the  time 
dependence  of  G^  through  the  contact  could  be  accounted  for.  In 
this  paper,  however,  we  shall  continue  to  ignore  this  variation 
as  the  essential  features  of  the  stress  variation  and  the  tractive 
force  behavior  are  illustrated  rather  well  without  it;  at  this 
point  including  it  would  seem  to  be  a somewhat  uncalled  for  re- 
finement. 

C.  The  Free  Volume.  While  the  computations 
reported  in  this  paper  are  carried  out  using  the  fictive  pressure 
concept  as  a means  of  characterizing  the  structure,  an  alternative, 
and  perhaps  more  familiar  approach,  makes  use  of  the  idea  of  free 
volume8  for  describing  the  structure.  The  dynamical  structure  para- 
meter £(t)  is  taken  to  be  the  instantaneous  free  volume 

vf(t)  » V(t)  - VC  (58) 


’ » n • • * « 


(58) 


21. 


where  V(t)  is  the  instantaneous  molar  volume  of  the  system  and  Vc 
is  the  "close- packed"  volume.  For  a system  in  equilibrium  at  a 
pressure  PQ,  the  time  dependence  of  the  structure  following  the 
imposition  of  a small  pressure  change  AP  at  time  ti  is  then  given 
by  an  equation  of  the  form  of  (40) ; specifically 

vf(t)  = vfQ  + (vfa8-vfo)  + (1-  -£=>  (vfoo  - vfo)  (l-4>  [ (t-ti)/T]  } . (59) 

In  this  equation  VfQ  and  Vf^  are  the  equilibrium  values  of  the  free 
volume  for  pressures  P0  and  PQ  + AP  respectively,  K*>/<  (which  plays 
the  role  of  Y)  is  the  ratio  of  the  instantaneous  (or  glassy)  and 
the  equilibrium  compressibilities.  The  generalization  of  this  to 
obtain  the  response  to  an  arbitrary  pressure  change  imposed  at  t^_ 

is  of  the  form  of  Eq. (50) : 

* 

t 

vf(t)  = vf0  - <«  /t]_  dt'  v(t')  (dp/dt') 

t t 

-(<»  “ <o^tl  dt'  v(t'}  (dP/dt’>  fl-'H/t'  dt"/ T (t " ) ] } , (60) 

' 

where  we  have  assumed  that  Vc  is  pressure  independent.  In  order 
to  use  this  we  must  relate  the  time  Varying  shear  and  structural 
relaxation  times,  ts  and  t respectively,  to  the  free  volume.  Gene- 

9 

rally  one  takes 

t or  ts  * exp  (A/vf) 

I 

where  A is  a constant.  Using  this  Eq.  (60)  provides  Vf(t)  which 
can  in  turn  be  used  to  give  the  time  dependence  of  rs.  This  last 
is  substituted  in  Eq.  (57)  to  give  the  effective  time  dependent 
viscosity  or  shear  stress. 
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III  DESCRIPTION  OF  THE  CALCULATION 

The  calculational  program  presented  in  the  preceeding  section 
has  been  carried  out  to  determine  which  of  the  viscoelastic  and 
material  properties  are  of  prime  importance  in  determining  the  time 
dependent  shear  stress  and  thus  the  traction.  We  consider  that  the 
pressure  variation  through  the  contact  zone  is  of  the  Hertzian  semi 
elliptical  form  [Eq. (12) ] : 

P(t)  = PHZ  [l-C2t/t0)2]1’  (04t«to)  . (61) 

For  the  maximum  Hertzian  pressure  PH2,  values  of  7.5  kbar,  15  kbar, 
and  30  kbar  have  been  considered.  Transit  times  t in  the  range 
10”^s  to  10~^s  (corresponding  to  rolling  speeds  larger  than  about 
80  cm/s  for  15  cm  diameter  steel  disks)  have  been  examined. 

To  describe  the  instantaneous  structure  of  the  lubricant  we 
have  used  the  concept  of  fictive  pressure  outlined  above.  The 
dependence  of  the  structural  and  shear  relaxation  times  on  the  in- 
stantaneous state  of  the  system  (described  by  P and  Pf)  is  assumed 
to  be  of  the  form 

(ts)  = (tsq)  exp  {a  [xP  + (l-x)Pf]>  (62) 

where  a is  the  pressure  viscosity  coefficient  (we  shall  assume  that 
G,,,  = ti/ts  is  pressure  independent)  . Observe  that  in  equilibrium 
(Pf  = P)  this  corresponds  to  the  familiar  n“exp(ap)  form  for  the 
pressure  variation  of  viscosity.10  The  "zero  pressure"  shear  re- 
laxation time  xso  is  given  in  terms  of  nQ  the  viscosity  at  the 


inlet  to  the  contact  zone  as 
tgo  = Hq/Goo 


(63) 
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and  the  ratio  of  structural  to  shear  relaxation  times  is  assumed  to 
be  a pressure  independent  characteristic  parameter  for  a given  sub- 
stance. Values  of  x generally  fall  in  the  range  0.4  to  0.6.^"^ 

The  time  evolution  of  the  fictive  pressure — given  by  Eq. (52) 
with  t]_=0 — and  thus  of  the  relaxation  times  is  governed  by  the 
structural  relaxation  function  (J>(t/-r).  The  shear  stress  dynamics, 
Eq.  C57)  , are  determined  by  the  shear  relaxation  function  <j>s(t/TS). 
Three  forms  of  these  functions,  which  we  distinguish  as  <j>  ^ (t/T)  , 

♦ ^(t/x)  and  <t>^(t/T),  were  considered  in  the  calculations.  The 
first  of  these  is  a -simple  exponential  decay 

<{>(1)  (t/T)  = exp  (-t/t),  (64) 

which  we  examine  in  the  spirit  of  a "zeroeth  order  approximation." 

A somewhat  more  realistic  form  is  the  so-called  fractional  exponen- 
tial decay  function 

<J>  (2)  (t/t)’  = exp  [-(t/T')B]  (0«B*1);  (65a) 

thf.'  average  relaxation  time  and  the  decay  constant  t'  are  related 
by 

x « dt<f> {2)  (t/x)  = (t'/B)  r (1/0)  (65b) 

where  T(  ) designates  the  gamma  function.  Experimentally  this 

form  is  found  to  do  a rather  good  job  of  describing  structural 

relaxation  data  in  a number  of  non-polymeric  substances  with  a value 

14-16 

of  B in  the  range  of  roughly  1/3  to  2/3.  In  polymeric  liquids 

the  short  time  part  of  the  response  is  also  fit  rather  well  by 
the  fractional  exponential  function  but  for  long  times  Ct>x'),  the 
decay  is  considerably  less  rapid;  that  is,  there  is  a long  time 
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"tail"  on  the  relaxation  function.  To  examine  the  behavior 

of  such  materials  we  consider  a third  form  of  the  relaxation 
functions: 

<f>(3)(t/-r)  = a exp  [-Ct/xJ5]  + (1-a)  exp  (-t/r")  (66a) 

where  0<a<l  and  t">>t'.  The  average  relaxation  time  for  this 
function  is  given  by 

x=  a(T'/S)rCl/S).  + (1-a)  t".  (66b) 

The  specific  functions  <f>  ^ and  <f>  ^ examined  are  characterized  by 
the  values  0 = Jj,  a = 0.8  and  t"/t'  = 200.  A comparison  of  the 

three  functional  forms  is  shown  in  Figure  6. 

- * . % 

Values  of  the  various  parameters  describing  the  viscoelastic 
behavior  were  chosen  in  ranges  typical  of  what  is  commonly  found 
in  lubricating  fluids.  In  addition  to  the  three  types  of  relaxation 
functions  we  take  the  following: 

Shear  Rigidity  Modulus  = Go  = 2 x 10^°  dynes/cm2; 

Structure  Coefficient  = x * 0.5 

Inlet  Viscosity  = n0  = 0.1P,  10P,  and  102P; 

Relaxation  Time  Ratio  = t/ts  = 1,3,  and  10 

Pressure  Viscosity  Coefficient  = a = 1 kbar“^,  2 kbar“l,  and 
4 kbar”^-  . 

The  actual  calculations  of  Pf(t)  and  y (t)  via  Eqs. (52)  and 
(57)  was  carried  out  by  dividing  the  time  range  t=0  to  t=t$  into 
N=200  subintervals  as  described  in  section  II. 4 above  (taking  a 
larger  number  of  subintervals  did  not  significantly  alter  the  re- 
sults) . The  intervals  were  selected  according  to  Eq.  (45)  such 
that  equal  magnitude  pressure  changes  AP  = ± 2PH2/N  occur  during 


25 


each;  the  pressure  during  the  k—  interval  (k<N/2) is  just  kAP 
(PQ  being  taken  as  zero) . The  relaxation  times  t and  ts  for  each 
interval  were  calculated  using  Eq.  (62)  in  which,  for  the  k—  in- 
terval, for  example,  P was  taken  as  kAP  and  Pf  was  assigned  its 
value  at  the  end  of  the  (k-1)— , i.e.,  the  preceeding,  interval. 
Thus,  for  instance,  during  the  first  interval  0<tst2,  the  relaxation 
times  are  related  to  their  "zero  pressure"  values  tq  and  tso  by 
the  factor  exp(ctxAP),  the  fictive  pressure  being  zero  (=PQ)  at  the 
start  of  this  interval.  Using  this,  Eqs.  (52)  and  (57)  can  be 
evaluated  straightforwardly  to  give  Pf(t)  and  u(t)  for  0^t<t2. 

For  the  next  interval,  the  relaxation  times  are  given  by  the  factor 
exp  { a [2xAP  + (1-x)  Pf(t2)]>  and  thus  Eqs.  (52)  and  (57)  provide 
Pf(t)  and  u (t)  for  t2<t^t3.  This  procedure  is  iterated  until  Pf(t) 
and  u(t)  are  mapped  out  over  the  entire  range  t=0  to  t=tQ.  It  is 
then  a straightforward  matter  to  integrate  y(t)  and  thus  obtain  the 
tractive  force. 

IV.  RESULTS  AND  DISCUSSION 

1.  General  Nature  of  the  Results 

To  illustrate  the  general  features  of  the  calculation  on 
which  our  attention  will  be  focused,  we  consider  the  results  of  a 
computation  performed  for  a maximum  Hertzian  pressure  PH2  35  15  kbar 
and  tQ  = 10~4s  (U  = 20  m/s) . The  characteristic  liquid  properties 


were  taken  as: 
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G,,,  * 2 x IQ10  dynes/cm2; 

nQ  = 10  P? 
x * **; 

a = 2 kbar”l; 
t/ts  = 3; 

= exp  C-t/f) 

♦ gCt/Tg)  = exp  C-t/Ts) 

The  computed  behavior  of  Pf(t)  and  y (t)  under  these  circumstances 
are  shown  in  Figure  7.  Consider  first  the  fictive  pressure  variation. 
Initially  the  viscosity  and  hence  the  relaxation  times  are  small  and 
the  system  is  able  to  equilibrate  in  response  to  the  applied  pres- 
sure changes.  As  the  pressure  becomes  larger  the  relaxation  time 
becomes  longer  so  that  the  system  is  no  longer  able  to  respond  to 
the  changing  pressure;  the  fictive  pressure  changes  relatively 
little  for  t £ 5 x 10"®s.  indicating  that  the  structure  of  the 
material  has  been  "frozen  in."  Observe  that,  from  Fig.  8,  this 
transition  from  viscous  liquid  to  amorphous  solid  behavior  sets  in 
when  t = 10“^ s,  that  is,  when  T is  on  the  order  of  the  time  scale 
for  the  occurrence  of  the  significant  pressure  changes.  Further 
increases  in  pressure  have  little  effect  on  the  materials ' struc- 
ture; for  most  of  the  time  that  is  is  in  the  contact  zone,  the 
lubricant  behaves  as  an  amorphous  solid  with  a structure  that  the 
equilibrium  material  would  have  at  a pressure  of  about  4 kbar. 
Eventually  for  t * 90  x 10  ^s,  the  pressure  has  fallen  to  a point 
where  the  relaxation  times  are  sufficiently  short  that  the  material 
can  again  respond  to  the  changing  pressure.  This  is  reflected  in 
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the  behavior  of  the  fictive  pressure  as  it  attempts  to  equilibrate 
at  the  value  of  the  applied  pressure. 

The  time  dependent  shear  stress  or  effective  viscosity  also 
reflects  this  kind  of  behavior.  Initially  (t  < a few  ys)  the  shear 
relaxation  time  is  quite  short  and  the  stress  is  relaxed  by  viscous 
flow.  However  when  t becomes  relatively  long  (t  > 10  ys)  the 
stress  increases  linearly  with  time  with  a slope  of  G*,;  this  is  just 
the  behavior  of  an  ideal  elastic  solid  for  which 

o (t)  = Goo  e(t)  = Goote;  y(t)  = G^t.  ’ (67) 

This  behavior  is  shown  by  the  broken  line  in  Fig.  7. 

From  Eqs.  (27)  and  (3)  the  traction  coefficient  can  be  written 
as  (for  small  slip) 


8R*  r I. 


Ct  = ty/W.t0]-G0  ~ 


l£  is  apparent  from  this  equation  that  a large  value  of  the 
shear  rigidity  Go,  leads  to  large  coefficients  of  traction.  It  is 
also  clear  that  (for  a given  film  thickness)  the  initial  slope  of 
the  traction  curve  is  determined  by  the  quantity  in  square  brackets . 
We  designate  this  as  x,  the  reduced  traction  slope: 

X = U/hG  tQ.  (69) 

For  the  situation  just  described,  integrating  the  y (t)  curve 
shown  in  Fig.  7 leads  to  X = 0.91;  the . initial  traction  curve  slope 
is  91%  of  what  would  be  obtained  for  an  elastic  solid  film  of  the 
same  thickness.  In  the  following  sections  we  consider  the  effect 
on  x of  varying  PHZ,  tQ,  nQ,  a,  t/ts  and  the  forms  of  4>  and  $>s.  In 


'*  ' »*  ,>v  + 


28. 


considering  the  effect  of  changing  PHZ,  tQ,  nQ  and  a we  should 
also  consider  the  variation  of  the  film  thickness  with  these 
quantities.  Empirically^ 

h « ct0-6  Cno/to)0'7PHZ0*45  (70) 

Consequently  we  shall  also  consider  the  quantity 

X*  S (hQ/h)  X (71) 

where  hQ  is  the  film  half-thickness  for  a * 2 kbar-!,  no  = 

PHZ  = kbar'  and  fc6  * 10_4s. 

2.  Effect  of  Varying -the  Rolling  Speed  and  the  Maximum 
Hertzian  Pressure 

The  time  variation  of  the  fictive  pressure  and 
shear  stress  were  computed  for  maximum  Hertzian  pressures  of  7.5  kbar 
for  several  values  of  the  transit  time  tQ.  The  reduced  traction 
slopes  X (as  well  as  the  values  of  X’)  were  obtained  from 

these.  The  lubricant  parameters  were  taken  as  a = 2 kbar-!, 
nQ  ■ 10  P and  f/Ts  58  3.  The  fractional  exponential  from  of  the 

g 

relaxation  functions,  e.g.,  $(tA)  = exp[-(t/r)  J,  with  8 = 0.5 
was  assumed.  The  results  are  presented  in  Figs.S  and  10  and  in 
Table  I.  The  general  behavior  is  quite  similar  to  what  was  noted 
above.  The  fictive  pressure  first  rises  rapidly  in  response  to  the 
impressed  pressure  change,  then  levels  off  when  the  structural  re- 
laxation time  becomes  long  and  ultimately  drops  precipitously  when 
the  pressure  is  reduced  as  the  lubricant  leaves  the  contact  zone. 

For  most  of  the  transit  time  the  structure  (the  fictive  pressure) 
is  unchanging  and  can  be  characterized  by  the  maximum  value  P^  max. 
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Observe  that  as  the  transit  time  is  decreased  (rolling  speed  in- 
creased) the  value  of  Pf  max  decreases  for  a given  value  of  PHZ . 

This  is  to  be  expected  since  for  small  transit  times  the  ability 
of  the  material  to  equilibrate  with  the  changing  pressure  is  in- 
hibited. 

Note  also  that  for  a given  value  of  tQ,  increasing  PHZ  causes 
a decrease  in  Pf  This  result,  surprising  at  first,  can 

be  understood  in  much  the  same  sort  of  way.  Increasing  PHZ  (for 
fixed  tQ)  also. increases  the  rate  at  which  the  pressure  is  changing 
and  consequently  reduces  the  ability  of  the  structure  to  remain 
in  equilibrium.  The'  effect  of  either  increasing  the  Hertzian  pres- 
sure or  decreasing  the  transit  time  is  to  cause  the  lubricant  to 
act  in  a solid-like  fashion  for  a larger  fraction  of  the  time  that 
it  spends  in  the  contact.  This  is  evident  in  the  behavior  of  the 
shear  stress  or  time  dependent  viscosity  u(t)  and  is  reflected  in 
the  variation  of  x with  both  tQ  and  PHZ . At  the  highest  pressure 
one  finds  essentially  elastic  solid  behavior  whereas  at  7.5  kbar, 
the  value  of  x varies  from  about  one-quarter  to  one-half  of  what 
would  be  found  in  a solid.  It  is  interesting  to  observe  that  if 
one  were  to  interpret  the  corresponding  traction  curve  slopes  as 
the  shear  modulus  G*  (effectively  putting  x * 1),  one  could  infer  as  much  as 
a four-fold  increase  in  Goo  as  PHZ  increases  from  7.5  kbar  to  30  kbar. 
This  is,  or  course,  entirely  spurious  as  we  have  taken  Go,  to  be 
constant  (2  x lO1^  dynes/cm2) . 


3.  Effect  of  Varying  the  Pressure-Viscosity  Coefficient 
To  examine  the  effect  on  traction  of  different  values  of  the 


pressure,  we  have  performed  the  calculations  for  the  situations 
outlined  below: 

External  Variables:  P„_  = 15  kbar 

iiZ 

tQ  = 100  pS  (U=20  m/s) ; 


Lubricant  Properties :G_  = 2 x 10  ^ dynes/cm^, 

no  = 10  P, 
t/ts  = 3, 

" <j>  = exp  [—  (t/ t ' ) * 5 1 » 

<f>s  - exp  [- (t/T s ' ) *5] , 

and 

a = 1 kbar"1,  2 kbar"1,  4 kbar  1. 

The  time  dependent  quantities  Pf(t)  and  \i  (t)  are  shown  in  Figure  11 
for  these  circumstances.  A summary  of  the  results  is  given  in  Table  II. 
From  these  it  is  apparent  that  increasing  a has  an  effect  quite  si- 
milar to  increasing  PH2 . This  is  not  surprising  as  a and  P appear  in 
the  calculation  as  a product  in  the  formula  for  the  pressure  de- 
pendence of  t and  ts;  consequently  their  influence  on  the  traction 
curves  is  symmetric. 

4.  Effect  of  Varying  the  Inlet  Viscosity 
For  the  situation  described  by  the  external  variables 
PH2  = 15  kbar, 
tQ  « 100  us  (U  = 20  m/s) 
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and  nQ  **  0.1P,  10  P and  10 3 P, 

the  results  are  shown  in  Figure  12  and  Table  III.  As  is  apparent, 
a higher  inlet  viscosity — corresponding  to  longer  relaxation  times 
at  the  entrance  to  the  contact — limits  the  ability  of  the  lubricant 
to  adjust  its  structure  to  the  pressure  changes  in  the  contact  and. 


tion  of  the  transit  time.  Observe,  from  the  X'  variation,  that  the 
major  effect  is  associated  with  changes  in  the  film  thickness. 

4.  Effect  of  Varying  the  Relaxation  Time  Ratio  t/ts 
For  the  same  external  conditions  as  were  considered  above  and 
for  the  liquid  properties  as  above  (with  no  = 10  P)  we  consider  the 
changes  that  result  from  allowing  t/ts  to  take  on  the  values  1,3  and 
10.  As  is  evident  from  the  results  given  in  Table  IV,  the  effects  on 
the  traction  parameters  x and  x'  are  quite  small.  The  decrease  in 


Pf/Inax  as  t/ts  is  increased  is  traceable  to  the  larger  values  of  t 
at  the  inlet  (ts  at  the  inlet  is  the  same  in  all  cases )•  consequently 
the  solid-like  behavior  sets  in  somewhat  earlier  for  larger  values 
of  t/ts. 

5.  Effect  of  Changing  the  Form  of  the  Relaxation  Functions 
To  examine  the  importance  of  the  detailed  shape  of  the  relaxa- 
tion functions  <J»(t/x)  and  <f>s(t/Ts)  in  determining  the  lubricant  be- 
havior in  an  EHD  contact,  we  have  examined  four  separate  cases: 


Case  I:  <J>(t/x)  ■ 

♦ gCtAs)  = <t>(1)  (t/xs)  ; 


Case  IX:  <j>(t/x)  = $C2)(t/x) 


4>s  ^/Ts)  = ♦ (2  ) Ct/r  ) ; 
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Case  III:  $(t/x)  = <f>(2)(t/x] 

*sCt/Ts)=  *C3)  (t/Ts); 

Case  IV:  <j»(t/x)  = $(3)Ct/x) 

4>s (t/xs)  = (f.C3)  Ct/xs)  . 

The  functions  <p  ^ and  <p  ^ are  defined  in  Eqs.  (64)  -(66) 

and  have  been  plotted  in  Fig.  6.  For  these  functions  we  have  taken 
0 = **,  a = 0.8  and  x"  = 200  x'.  For  each  of  the  four  cases,  we  com- 
pute Pf(t),  y(t),  x and  X'  f°r  the  following  conditions: 

External  Variables:  PHZ  = 7.5  kbar, 

tQ  = 100  ys  (U  = 20  m/s)  ; 

Lubricant  Properties:  G^  = 2 x 10^-°  dynes/cm2 

Ho  = 10  p 
a * 2 kbar“l 

t/ts  = 3 

The  results  are  shown  in  Figure  13  and  Table  V. 

We  consider  first  the  behavior  of  the  fictive  pressure.  It  is 
clear  from  the  results  that  for  relaxation  functions  characterized  by 
the  same  average  relaxation  time  I (equal  viscosities  here  since  G* 
and  x/x_  are  the  same) , there  is  a rather  marked  difference  in  the 
behavior  of  the  fictive  pressure.  This  stems  primarily  from  the  en- 
hancement of  the  short  time  response  materials  characterized  by  a struc- 
tural relaxation  function  of  the  ^ type  as  compared  with  those  cha- 
racterized by  <p  ^2^  (which,  in  the  same  way  allows  for  more  short— 
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time  response  than  The  effect  of  this  short  time  response  ca- 

pability can  be  understood  by  thinking  of  the  relaxation  functions 
as  superpositions,of  exponential  decays  of  the  form 

E g exp (-t/T . ) , (72) 

i 1 1 

where  the  g.^  are  weighting  factors.  Suppose  for  definiteness  we  take 
x (the  average  relaxation  time)  = tq,  i.e.  the  i=0  term  in  the  sum, 
and  let  i run  from-M  to  +M,  with  x^  < x^+ 

For  the  simple  exponential  decay  function  only  gQ  ^ 0 and  we  may 
think  of  the  structure  as  responding  on  a time  scale  x0.  When  the 
pressure  changes  occur  on  a scale  long  compared  with  xQ  the  structure 
is  able  to  remain  in  equilibrium.  However  when  t increases  (because 
of  the  pressure  increase)  to  a value  comparable  with  the  time  scale  of 
pressure  changes  (for  the  case  here  when  xQ  = I0”5s) , the  material  is 
no  longer  able  to  adjust  its  structure  and  consequently  departs 
from  P. 

For  the  functions  <j>  ^ and  <j>  however,  all  the  g^_  are  non-zero; 
in  particular  there  are  short  relaxation  time  (x^<x)  contributions  to 
the  structural  relaxation  function — the  i<0  terms  in  (72) . Conse- 
quently even  when  To=i0_^s/  there  are  still  structural  relaxation 
processes  that  are  able  to  respond  and  as  a result  there  is  a con- 
tinuing (albeit  incomplete)  evolution  of  the  structure.  This  is  par- 
ticularly evident  in  the  fictive  pressure  curves  H, III  and  17 in  Fig. 13a. 

The  influence  of  the  short  relaxation  time  contributions  on  the 
shear  stress  (effective  viscosity)  behavior  is  somewhat  different. 

The  specific  nature  of  the  difference  can  be  understood  most  easily 
by  comparing  the  curves  II  and  III  in  Figure  13b.  -For  these,  the 
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structural  relaxation  Cbut  not  the  shear  relaxation)  dynamics  are 
the  same.  As  is  evident  the  shear  stress  is  significantly  less  (lead- 
ing to  smaller  traction  slopes)  when  the  short  time  shear  relaxation 
processes  are  enhanced.  Simply  stated  this  means  that  even  when 
those  contributions  to  the  stress  described  by  long  relaxation  times 
are  kinetically  arrested  and  are  thus  behaving  in  a solid-like 
manner  (again  we  are  thinking  of  representing  <j>s  as  a weighted  sum 
of  exponential  decays) , there  are  still  stress  components  with  short 
relaxation  times  present  that  can  be  relaxed  in  a liquid-like  fashion. 

This  partial  stress  -relaxation  produces  the  lower  stress  values  and 

(3) 

smaller  traction  slopes  for  the  case  <j>s  = 4> 


V.  SUMMARY 

The  influence  of  various  liquid  lubricant  properties  on  the  slope 
of  the  traction  curves  in  the  low  slip/roll  region  has  been  evalua- 
ted. This  slope,  given  by  an  expression  of  the  form 
(8R/TThE)  X G., 

is  quite  sensitive  to  certain  of  these  properties  and  relatively  in- 
sensitive to  others.  The  effects  of  each  are  summarized  below. 

1.  Hertzian  Pressure:  Increasing  PHZ  causes  x to  approach  unity 

and  (for  a given  rolling  speed)  causes  the  film  thickness  h to  de- 
crease; the  traction  slope  therefore  increases.  If  the  transit  time 
tQ  is  fixed,  increasing  PHZ  causes  b to  increase.  Since  this  in- 
creases U,  h is  increased  and  this  change  opposes  the  change  in  x* 

2.  Rolling  Speed:  Increasing  U causes  x to  approach  unity  and  pro- 
duces an  increase  in  h.  For  relatively  low  values  of  PgZ  the  ^onner 
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effect  dominates  and  the  traction  slope  increases;  however,  when  P 

HZ 

is  sufficiently  high  that  the  latter  effect  assumes  the  major 

role  and  causes  a reduction  in  the  traction  slope  with  increasing  U. 

3.  Shear  Rididity:  As  is  evident  an  increase  in  will  cause  an 
increase  in  the  traction  slope.  This  is  the  single  most  important 
factor  although  the  practical  range  of  Go,  values  that  are  available 
is  rather  limited. 

4.  Pressure-Viscosity  Coefficient:  Increasing  a causes  x to  increase 

toward  unity  but  also  produces  an  increase  in  h.  At  lower  values  of 
PHZ'  t*ie  increase  in  x predominates  and  this  the  traction  slope  in- 
creases; when  PH2'is-high,  xal  and  the  increase  in  h causes  the 
slope  to  decrease. 

5.  Inlet  Viscosity:  Increasing  nQ  produces  essentially  the  same 

changes  as  increasing  a.  However,  since  rather  large  ranges  of  n0 
values  are  available,  the  film  thickness  increase  is  of  greater  im- 
portance. 

6.  Relaxation  Time  Ratio  x/xs:  Increasing  x/xs  causes  a decrease  in 

x and  in  the  traction  slope.  For  the  ranges  of  x/xs  generally 

encountered  this  effect  is  quite  small. 

7.  Form  of  the  Structural  Relaxation  Function:  Increasing  the  short 
time  response  region  of  the  function  <J>(t/x)  leads  to  small  increases 
in  x and  thus  in  the  slope.  If  G,*  had  been  allowed  to  vary  through 
the  contact  region  (rather  than  being  held  constant) , this  increase 
would  have  been  more  pronounced  as  higher  fictive  pressure  values  and 
thus  higher  shear  moduli  are  achieved  when  a short-time  portion  of 

<$  (t/ t ) is  present. 


36. 

8.  Form  of  the  Shear  Relaxation  Function:  Increasing  the  short 
relaxation  time  poriton  of  $s(t/Ts)  leads  to  a decrease  in  x and 
therefore  to  a decrease  in  the  traction  slope. 

The  liquid  properties  of  chief  significance  in  producing  a 
large  reduced  traction  slope  x are  therefore  a large  shear  rigidity, 
relatively  large  values  of  the  inlet  viscosity  and  pressure-viscosity 
coefficient,  a structural  relaxation  function  with  an  enhanced  short 
time  response  nature  and  a shear  relaxation  function  without  short 
time  response  characteristics . 

The  influence  of  a pressure  and  structure  dependent  shear 
rigidity  and  the  possibility  of  extending  the  calculation  into  the 
non-linear  region  (which  is  of  greater  practical  significance)  of 
large  AU/U  values  are  presently  under  investigation. 
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TABLE  I.  The  effect  of  changing  the  rolling  speed  and  the  maximum 
Hertzian  pressure  . 


tQ  Cys) 

U (itv/s) 

p.  (kbar) 

rf ,max 

X 

X* 

PHZ  = 

7.5  kbar 

10 

100 

4.33 

0.523 

0.143 

40 

25 

5.04 

0.420 

0.320 

50 

20 

5.15 

0.402 

0.338 

66.7 

15 

5.29 

0.381 

0.392 

100 

• * 

10 

5.48 

0.349 

0.477 

200 

5 

5.80 

0.297 

0.659 

400 

— 2.-5 

6.10 

0.245 

0.883 

PHZ  = 

15  kbar 

10 

200 

3.56 

0.900 

0.180 

80 

25 

4.75 

0.851 

0.728 

100 

20 

4.87 

0.845 

0.845 

133 

15 

5.03 

0.835 

1.019 

200 

' 10 

5.26 

0.825 

1.340 

400 

5 

5.65 

0.805 

2.124 

PHZ  = 

30  kbar 

10 

400 

2.74 

...996 

0.146 

100 

40 

4.04 

.984 

0.720 

(a)  Liquid  parameters 

taken  as  GM  = 2 

x 10^®  dynes/cm2 

t/  * 3,  u **  2 kbar”^. 
Eq.  (65)  with  8 = . 5 . 


Relaxation  functions  are  given  by 


TABLE  III.  Effect  of  changing  the  Inlet  Viscosity. 


(a) 


lo  w 

pf ,max 
(kbar) 

X 

x' 

0.1 

7.42 

.694 

17.4 

10 

4.87 

.845 

0.845 

103 

2.26 

.942 

0.038 

= 2 kbar-^ , 

nQ  variable; 

other  parameters  as 

in  Table 




TABLE  IV.  Effect  of  changing  the  relaxation  time  ration  t/ts.  ^ 


T/Ts 

P- 

f ,max 
(kbar) 

X 

x' 

1 

5.48 

.868 

.868 

3 

4.87 

.845 

.845 

10 

4.19 

.812 

.812 

= 10  P, 

t/t  variable; 

S 

• 

other  parameters 

as  in  Table 

IL 


TABLE 


V.  The  effect  of  changing  the  form  of  the  relaxation 
function  $ and  $s. 


f ,max 


X‘ 


p and  $ 

”s 

P - P(1)  (t/T) 

,g  » (j,  (1)  (t/Tg  ) 

♦ - <P(2)  (t/T) 

>S  ■ ^ ^ (t/Tg) 

<P  * <P  (t/x) 

»g  - (t/Tg) 

<P  ■ <j>  (t/x) 

»s  * (3)  (t/xs) 


Ikbar) 

5.14  .444 


5.48  .349 


5.48  .194 


6.18  .225 


.607  Case  I 


.477  Case  II 


.265  Case  III 


.307  Case  IV 


External  variables:  PHZ  * 7.5  kbar,  tQ  = 10-4s?  liquid  properties: 
n0  * 10  Pf  o = 2 kbar-1,  x/xg  = 3. 

lb)  n»  #2)  (3) 

i p , <p  and  <f>  are  defined  in  Eqs.  (64)-(66);  we  have  taken 
0 * .5,  a = .8  and  x"  = 200  x'. 


i > #V 
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Figure  1.  Twin-disk  apparatus  for  traction  measurements. 

Figure  2.  EHD  contact  zone  geometry. 

Figure  3.  Typical  traction  curves  measured  on  a twin-disk  apparatus 
for  two  maximum  Hertzian  pressures  [Data  from  K.L. Johnson 
and  R.  Cameron,  Proc.  Inst.  Mech.  Engrs.  (1967-68) , 

182  Pt.  1,  307.] 

Figure  4.  (a)  A presssure  step  AP  applied  at  time  ti.  (b)  The 

• response  of  the  structure  to  that  pressure  step. 

Figure  5.  The  "series-of-steps"  approximation  (solid  curve)  to  a 
smooth  pressure  variation  (dash-dot  curve) . 

Figure  6.  The  relaxation  function  $(1*,  <^2)  and$(3),  defined  in 
Eqs.  (64)  - (66) , plotted  versus  reduced  time. 

Figure  7.  Response  of  a system  under  the  conditions:  Phz  = 15  kbar, 
tQ  - 100us,  rio  “ 10P,  a = 2 kbar”1,  t/ts  = 3,  and 
$(t/t)  = exp  (-t/t) . (a)  The  pressure  P (broken  curve) 

and  fictive  pressure  Pf  (solid  curve)  versus  time,  (b)  The 
effective  viscosity  or  reduced  shear  stress  versus  time 
(solid  curve) . The  behavior  of  an  ideal  elastic  solid 
is  shown  (broken  curve)  for  comparison.  The  arrow  at 
the  left  indicates  the  average  effective  viscosity  IT. 

Figure  8.  The  time  variation  of  the  structural  and  shear  relaxa- 
tion times  (t  and  ts  respectively)  for  the  same  con- 
ditions as  listed  in  Fig.  7. 

Figure  9.  The  effect  on  the  structural  and  shear  response  of  vary- 
ing the  transit  time  tQ  (tg  values  shown  at  the  right 
of  the  graphs) . (a)  Pressure  (broken  curve)  and  fictive 

pressure  (solid  curves)  versus  t/tQ.  (b)  Normalized 
effective  viscosity  or  shear  stress  versus  t/to  (solid 
curves) ; elastic  solid  behavior  (broken  curve)  is  shown 
for  comparison.  The  arrows  at  the  left  indicate  the 
average  effective  viscosity.  For  these  curves  Phz=7.5  kbar 
Ho  * 10P,  a * 2 kbar”1,  t/ts  ■ 3 and  <j>  and  <|>s  are  given 
by  Eq.  (65) . 

Figure  10.  The  effect  on  the  structural  and  shear  response  of 
varying  the  peak  Hertzian  pressure  %z  (Phz  values 
shown  at  the  right  of  the  graphs) . (a)  Fictive  pressure 

versus  time,  (b)  Effective  viscosity  versus  time  (solid 
curves) ; elastic  solid  behavior  (broken  curve)  is  shown 
for  comparison.  The  arrows  at  the  left  indicate  the 
average  effective  viscosity.  For  these  curves  the  con- 
ditions are  as  in  Fig.  9,  except  that  to  * 100  ys  and 


HZ 


takes  on  the  values  shown. 


Figure  11.  The  effect  on  the  structural  and  shear  response  of 

varying  the  pressure-viscosity  coefficient  a (a  values 
shown  at  the  right  of  the  graphs) . (a)  Fictive  pressure 

versus  time,  (b)  Effective  viscosity  versus  time  (solid 
curves) ; elastic  solid  behavior  (broken  curve)  is 
shown  for  comparison.  The  arrows  at  the  left  indicate 
the  average  effective  viscosity.  For  these  curves  the 
conditions  are  as  in  Fig.  10  except  that  PHz  = 15  kbar 
and  a takes  on  the  values  shown. 

Figure  12.  The  effect  on  the  structural  and  shear  response  of 

varying  the  inlet  viscosity  n0  (no  values  shown  at  the 
right  of  the  graphs).  (a)  Fictive  pressure  versus  time, 
(b)  Effective  viscosity  versus  time  (solid  curves) ; 
elastic  solid  behavior  (broken  curve)  is  shown  for  com- 
parison. The  arrows  at  the  left  indicate  the  average 
effective  viscosity.  For  these  curves  the  conditions 
are  as  in  Fig.  11,  except  that  a - 2 kbar-1  and  no 
takes  on  the  values  shown. 

Figure  13.  The  effect  of  different  structural  and  shear  relaxa- 
tion behavior  on  the  response  of  the  system.  The  four 
cases  shown  are  tabulated  at  the  beginning  of  Section 
IV-5.  (a) Fictive  pressure  versus  time.  (b)  Effective 

viscosity  versus  time  (solid  curves) ; elastic  solid 
behavior  (broken  curve)  is  shown  for  comparison.  The 
arrows  at  the  left  indicate  the  average  effective  vis- 
cosity. For  these  curves  the 'conditions  are  as  in 
Fig.  12  except  that  Prz.  » 7.5  kbar,  nQ  = 10P  and 
<|>  and  <j>s  take  on  the  forms  indicated. 
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